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I. INTRODUCTORY REMARKS 

The very possibility of applying the modern meth- 
ods of the classical theory of non-linear oscillations to 
quantum mechanics is based upon the representation of 
the non-stationary Schrodinger's equation as a classical 
Hamiltonian system. From this perspective it is quite 
natural to construct a special asymptotic perturbation 
theory which utilises the advantages of the Hamiltonian 
formalism, that is to apply canonical transformations. 
Special and sufficiently efficient approaches 0, and 
PI were developed by mathematicians for canonical sys- 
tems. However, the generality of these approaches makes 
them very cumbersome whereas the first two non-trivial 
approximations are ordinarily sufficient for practical ap- 
plication. 

Traditionally, mathematicians use the methods of 
spectral analysis of operators for constructing perturba- 
tion theory in non-relativistic quantum mechanics, see 
0,0. However, it is necessary to take into account that, 
in practice, physicists do not distinguish between the con- 
cepts of self-adjoint and symmetric operators. This gives 
rise to two unpleasant things: firstly, the domain of def- 
inition of the operator in Hilbert space remains unclear 
which does not allow one to apply the methods of spec- 
tral theory, and secondly, the domain of definition of the 
operator includes all functions for which analytical op- 
erations are meaningful regardless of the fact whether 
these functions (and the result of applying an operator 
to them) belong to Hilbert space. A rigorous consider- 
ation of the latter case requires the introduction of an 
equipped Hilbert space, see 0- 

For this reason it is not possible to prove even the 
conditions for applicability of the regular perturbation 
theory developed by Kato-Relih , U providing us with 
the criterion for the Rayleigh-Schrodinger formal series 
to have a non-zero radius of convergence. 

It is well known that orthodox perturbation theory is 
not applicable to many cases since the corresponding se- 
ries diverge. The asymptotic character of the series used 
in perturbation theory was first proved by Titchmarsh 
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[9j . A rigorous proof of the divergence of this series for 
an anharmonic oscillator {V ~ x'^) is given by Bender 
and Woo [13. 

In addition to the associated complexity, the above- 
mentioned methods possess another shortcoming, namely 
they do not allow one to obtain the wave function which 
plays an important role in the investigation of physical 
systems. 

In the present chapter we suggest another approach 
which is based upon the representation of the non- 
stationary Schrodinger's equation as a classical Hamilto- 
nian system. This representation enables one to make 
use of the powerful, modern, rigorously substantiated 
methods of the classical theory of non-linear oscillations 
(asymptotic perturbation theory) and indicate simple 
conditions for justifying the applicability of the results 
obtained. 

An important part for the transition from the hypothe- 
ses of Planck and Einstein to quantum mechanics was 
played by the adiabatic hypothesis by Ehrenfest. Born 
and Fock showed in 1928, that Ehrenfest's hypoth- 
esis is a consequence of the postulates of quantum me- 
chanics. A rigorous mathematical proof of the adiabatic 
theory was given by Kato in 1949, 0. Later on, the 
adiabatic Landau-Dykhne approximation, 01' E3i 
was built on the analogy between the adiabatic and quasi- 
classical approximations. 

The Born- Fock adiabatic approximation is actually not 
an approximation since all of the terms of the adiabatic 
Born- Fock series have the same order of smallness, (ISj, 
|16|. which, in turn, does not allow us to construct a 
post-adiabatic approximation. The Born-Fock condition, 
which implies real-valued wave functions, does not allow 
us to use this approximation in problems involving mag- 
netic fields. 

The results of the Landau-Dykhne adiabatic approxi- 
mation relate to the results of non-stationary perturba- 
tion theory only approximately. In addition to this, both 
approximations yield an incorrect factor in front of the 
exponential function, see p^ . 

The problem of the time interval, within which the 
difference between the approximated and exact solutions 
is small, plays an important part in the non-stationary 
case. In the above works, this problem is not discussed 
at all. 

In the present chapter, the adiabatic and the post- 
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adiabatic theories, as well as the adiabatic perturbation 
theory and the post-adiabatic approximation, are con- 
structed by means of the method of canonical averaging 
(phase perturbation theory). The basic assumptions of 
the Born-Fock theory are not satisfied. The obtained ap- 
proximations are compared with exact solutions of the 
non-stationary Schrodinger's equation for a harmonic os- 
cillator in a homogeneous time-dependent field and with 
the approximations obtained by traditional formulae |l3j | . 
One can see from this comparison that the standard 
non-stationary approximations [l3i | are valid only within 
a non-dimensional time intervals t ~ 1, whereas the 
approximations of the present chapter are valid within 
asymptotically longer time intervals t ^ 1/e. 



II. STATIONARY SCHRODINGER'S 
EQUATION AS A CLASSICAL HAMILTONIAN 
SYSTEM 

In this section the classical canonical perturbation the- 
ory is applied to constructing asymptotic solutions of 
Schrodinger's equation with a discrete spectrum. The 
main subject of analysis of the non-relativistic quantum 
theory is Schrodinger's equation, [T^ . 



ih 



dt 



(1) 



where = — 1, h — 1.054 10"^'* Js denotes Planck's 
constant, q = (91,92, ■•■,9n) denotes a point of the con- 
figuration space of the corresponding classical system, t 
is time and Vl/(<7, denotes a complex- valued function 
with integrable square of the absolute value. Further, 
H denotes a self-adjoint (symmetric) operator in Hilbert 
space, which in Cartesian coordinates, in Schrodinger's 
representation for one particle, has the form, |13| . 



H^T + V^~—A + V{x,y,z), 
2m 



(2) 



where f and V designate operators of the kinetic and 
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potential energies, respectively, A = -— ^ -I- — ^ -t- tt-^, 

ox'^ oy^ oz'^ 

and m is mass of the particle. Schrodinger's equation 
(0 is subject to an initial condition ^'(9,0) — ^'o(<z) and 
some boundary conditions. 

The case studied in the framework of perturbation the- 
ory appears when operator H can be cast as the sum 



H^Ha+eV, < e << 1, 



(3) 



of two self-adjoint operators, the corresponding problem 
(|H) for operator Hq being assumed to have an exact solu- 
tion and the second operator (perturbation) being small 
in some sense, [T^j . 

The majority of the physically interesting problems 
turn out to be mathematically incorrect, since the per- 
turbation operators are usually not boun-ded and not 



even self-adjoint. The latter is related to the fact that 
physicists never distinguish between the concepts of self- 
adjoint and symmetric operators. This leads to the oper- 
ator space in Hilbert space being unclear, which in turn 
does not allow one to apply the methods of spectral the- 
ory 0) Hi' For this reason, it is difficult to indicate 
the conditions for the applicability of perturbation the- 
ory and estimate the discrepancy between the exact and 
an approximate solution for practical problems. 

However, it is possible to reduce Schrodinger's equa- 
tion to a form of classical Hamiltonian system, which is 
well-developed in non-linear mechanics. This allows one 
to apply the methods of classical dynamics which are rig- 
orously substantiated and simpler from the perspective 
of application. 

Let us consider eq. with Schrodinger's operators 
(O, i.e. the problem 



in — — — = (^0 



dt 
*(9,0) 



eV)^{q,t), 



*o(g), 



(4) 



where operator Hq does not depend on time and e is a 
formal small parameter. The question of choosing the 
small parameters is discussed below. 

Along with the problem we consider, the generating 
approximation, which is obtained from eq. Q at e = 



ih 



5*°(q,i) 
dt 



^Ho^°{q,t), ^"{q,0)^^o{q)- (5) 



Assuming the spectrum to be discrete, we can apply 
Fourier's method and set the general solution of prob- 
lem Q in the form 

00 

vl/°(g,t) = ^c°^°(<7)exp(-zc.°0, 

n=0 

cl = [Mq)i^niq)dq, u;l^E°Jh, (6) 



where ipniq) ^^'^ denote respectively the eigenfunc- 
tions and eigenvalues of the following problem 



Ho^liq)^Ey^M 



(7) 



and an asterisk denotes the complex conjugate. 

The self-adjoint character of operator Hq means that 
for '^{q,t) S the following expansion is valid 



^{q, t) = J2 Cnmliq) exp(-zwOO 



(8) 



Inserting expansion ((SJ into eq. Q yields the following 
equations for the coefficients of the expansion c„(t), c* (i) 



Cn{t) = -ie ^ u„mCm(i)exp(-ia;^„i), 

m=0 
oc 

c*n{t) = ie^Vmncl,{t)eyip{iu:l,J), (9) 



m=0 
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vmuit) ^j_J V^*(g)F(g,t)^°(g)dq, (10) 

where a dot implies time derivative and = oj^^ — wJJ . 

The system of equations (jSjl is Hamihonian (in the 
classical sense) with the following Hamilton function 



eHi{c,c*,t) = -ie ^ i;„™c* c„ exp(iw°„t) (11) 

n,m— 

and describes a classical distributed system with an infi- 
nite number of internal resonances. The system is Hamil- 
tonian as matrix Vnm is Hermitian, that is the perturba- 
tion operator is self-adjoint. By separating the principal 
resonance uj^ — uj^ , we can cast Hamilton's function Hll|) 
as follows 



eHi{c, c*,t) 



oo 



-ie ^ 'w„mC* Cm exp(iu;°^t) (12) 



n,m— 



where a prime denotes the sum without the term with 
n = m. Transformation of variables c„,c* to the real- 
valued " action-angle" variables /„ , tpn by means of the 
formulae 



'l^exp{iipn),'fpn = - arctan " (13) 



«(Cn +c;^) 



we can set system in the form 



I71 = £2 ^ Wlnlrr 



m— 



X Im{u„„exp[-i(V'm - V'« + wJJ^n*)]} 



(14) 



m— 

X Re{w„™exp[-i(V'm - i>n + Wm„i)]} 
and Hamilton's function p2ll as follows 

00 

m.,n=0 (15) 

X ea::p[-i(V'm - "iAn + ^^m„i)] 

Systems (O and H14|l do not contain Planck's constant 
explicitly and are the classical Hamiltonian systems with 
an infinite number of internal resonances. Estimates of 
the norm of discrepancy between the exact and approxi- 
mate solutions as well as the conditions for applicability 
of the averaging method for these systems are given by 



the Los theorem, see [13 and HI, which is a generali- 
sation of Bogolyubov theorem for the case of an infinite- 
dimensional coordinate Hilbert space. 

The canonical form of systems (jHl and (|14|l allows us 
to consider the evolutionary equations by operating only 
with Hamilton's functions H12I) and IjlSfl . i.e. by calculat- 
ing an averaged Hamilton's function. For example, for 
eqs. (jllll and (|15|l the second approximation ff^^^ for 
the averaged Hamilton's function is constructed with the 
help of the following formulae 



if (2) (c, c* ) = eHi (c, c*) -I- £2^2 (c, c*), 



dc 



dc 



(16) 



where c„ and c* denote the evolutionary components of 
variables c„ and c*. In the latter equation the following 
notation is used 



to+T 



(/) = hm - 



fic, c*,t)dt, 



/(c, c*,t) = /(c, c*,t)-(/), 
{/} = / />, c*,t)dt, 



(17) 



the arbitrary function of slow variables c and c* being set 
to zero while estimating the last integral. 

Let us notice that H is an integral of the averaged 
equations of motion, i.e. an adiabatic invariant pol |. |2l| . 

The first approximation Cn ^ to expansion coefficient c„ 

is given by the formula cl^^ = c„, where c„ satisfies the 
following equation 



dc* 



(2) 

The second approximation c„ to expansion coefficient 
c„ is given by 



(18) 



f:(2) 



c„ + e- 



d{Hi} 
del ' 



(19) 



where the second approximation to evolutionary compo- 
nent c„ is obtained from the equation 



dct dct 



(20) 



III. GENERAL PROPERTIES OF THE 
CANONICAL FORM OF SCHRODINGER'S 
EQUATION 

Representation of Schrodinger's equation in canonical 
form or H14|) allows us to draw a number of conclusions 
without performing any calculations. 
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1. The original formulation of the problem of perturba- 
tion theory for the non-stationary Schrodinger's equation 
results in formulae enabling us to study all cases: sta- 
tionary (non-degenerate and degenerate), non-stationary, 
resonant, adiabatic etc. 

2. In the stationary case, for which matrix elements 
Vmn are independent of time, in the absence of degener- 
acy and the internal resonances 7^ O (e) , m ^ n) 
except for the selected principal resonance {m = n), sys- 
tems © or (I14II admit direct averaging. The equation 
for the evolutionary components of phase ipn has the 
form ipn = £Vnn (the mean values of the sums in eqs. 
© or l|T^ are zero) and is easily integrated, to give 
ipn = £Vnnt + 4'nO- In turn, this means that the first 
correction to the energy appears due to the principal in- 
ternal resonance. This conclusion remains valid for any 
order of the perturbation theory since the corrections to 
the eigenvalues of the unperturbed operator Hq are ba- 
sically determined by the principal internal resonance. 

On the other hand, it is evident that the coefficients 
of the diagonal terms in the averaged Hamilton function 
are the corrections to the energy. 

Let us notice that direct averaging without introduc- 
ing the resonant terms would lead to the appearance of 
divergent terms which are proportional to 1/ (£"„ — £"„) 
(to — > n). Thus, separation of the resonant terms elimi- 
nates the divergent terms in the perturbation theory se- 
ries. 

3. It is obvious that the first-order contribution to 
the averaged Hamilton function or to the above equa- 
tions from the non-resonant sums is zero for a stationary 
perturbation and non-zero for a non-stationary perturba- 
tion. This is the only difference between the stationary 
and non-stationary cases. In both cases, calculations are 
carried out by means of formula Ijltill . 

4. Contributions from these sums are not zero if the 
perturbation is stationary and the unperturbed levels 
contain close levels for which uj^^ — 0{e). Thus, the 
problem of close levels should be solved only in a non- 
stationary form. It is evident that the problem for de- 
generate levels is a particular case of the previous one 
for which, along with the principal quantum number n, 
there is a multiindex a characterising the unperturbed 
eigenvalue, so that the relationship ujna,nf3 = O (e) holds. 
All these cases are manifestations of the resonances ad- 
ditional to the principal internal resonance. 

5. Generally speaking, a general analysis is not ap- 
plicable for time-dependent perturbations because it is 
necessary to know the spectrum of the perturbation so 
as to take correct account of the possible resonances. 
Only general analysis of periodic (in particular, single- 
frequency) perturbation is feasible. 

6. For the sake of simplicity, let us consider the 
case of single-frequency perturbation {vmn ~ exp(±jwt)). 
Clearly, the problem in this case is reduced to the station- 
ary one, for which the expansion is performed for wave 
functions of the stationary states with the new frequen- 
cies LOmn = — W„ SUCh that LOm = l^m ~ ^ 1^^ ^ri = 



-I- w/2 (for V ~ exp(iwf)) and = -f a;/2, w„ = 
— Li;/2 ( for w ~ exp(— zwt)). 

Thus, the corresponding quantities = ^m,n are 

energies, that is the energy of the system: unperturbed 
system (atom) and a field. Their interaction is absent, 
the original interaction being included into the definition 
of the quasi-energy. This gives rise to the concept of the 
system dressed by a field (dressed atom, HJ)- 

Traditionally, this conclusion is obtained in a rather 
sophisticated way by means of Floquet theorem and the 
conclusion on the level splitting is obtained only in the 
resonant approximation, see p^. pol | . Let us notice that 
the results of this point are valid for any value of param- 
eter e, that is regardless of the perturbation theory. 

7. The exact eigenfrequencies of the perturbed sys- 
tem are given by the relationship — ilP^-V ^n- Thus, 
the second equations in (|14|l determine corrections to 
the eigenfrequencies caused by perturbation. In prin- 
ciple, these corrections can be removed by a canonical 
transformation of the phase which can serve as a dis- 
tinctive procedure of renormalisation which allows one to 
remove the secular terms from the series of the perturba- 
tion theory. It is obvious that in the non-stationary case 
(even for w„„ = 0) there exists a non-zero contribution of 
the first order stemming from the sum. The presence of 
this contribution is not possible using orthodox perturba- 
tion theory. The coherent interaction with the external 
field is realised under the condition of constant phase 
difference (the condition of equality of the original fre- 
quencies is only a necessary condition) and has the form 

¥'2 — </5l — = + "02 ^ ^1 + "01 — ^21 ~0J + '>p2—1pl — 0, 

where uj denotes the frequency of the external field sat- 
isfying the condition that cj^i — w = O (e). 

One usually uses the condition that the transition fre- 
quency is close to that of the external field. It follows 
from the form of the equations for tpn that the type of 
problem for resonant interaction coincides with the type 
considered in point 4 above. 

The exact frequencies (j;„n the system are 
anisochronous which is a characteristic of the non- 
linear classical system and leads to bounded solutions at 
resonance even in the case of no damping, in spite of the 
linearity of Schrodinger's equation. 

Thus, Schrodinger's equation is equivalent to some 
classical non-linear distributed system whereas represen- 
tations (PI or H14I) are expansions in terms of normal 
forms of the unperturbed system. 

To some extent, it is the picture to which Schrodinger 
tended and which is most close to the classical one. 
"There is no need to explain that the representation of 
the energy transformation from one oscillatory form to 
another under a quantum energy transition is much more 
satisfactory than that of electron jump"^, 23]. 

8. It is evident that no specific " quantum- mechanical" 
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properties of the generating operator Hq, but being self- 
adjoint, is used. Nevertheless, this technique allows us 
to analyse other self-adjoint problems of mathematical 
physics with a discrete spectrum. Taking into account 
the particular structure of the generating operator Hq 
we can construct a more efficient theory 24] . 

An attempt to apply the averaging method to quantum 
mechanics was undertaken in |2a|. However, as follows 
from the above, the absence of resonances (the main as- 
sumption of the authors) is not correct. Besides, despite 
the title of this article, the authors did not succeed in a 
proof of the theorem of convergence of the constructed 
perturbation theory. 

9. All formulae remain valid for the case of adiabatic 
perturbation, i.e. under the additional dependence of 
perturbation on the slow time t — et {eV — eV {r, t,T)). 
In this case, equations for the evolutionary components 
become non-stationary and require more sophisticated 
integration methods. 

10. Small parameter e is introduced in systems © and 
(|14|l in a formal way. Generally speaking, the question 
of a rigorous introduction into equations should be con- 
sidered individually for each particular problem. Let us 
point out some general ideas. 

Let us introduce some characteristic values [E] and [V] 
for the eigenvalues E and matrix elements Vmn respec- 
tively. Then [E] /H = luq and [V] /h = Slo can be referred 
to as a characteristic eigenfrequency and the generalised 
Rabi frequency (for a dipole interaction V ~ (IEq and 
V/h ^ n is called the Rabi frequency). Further study 
depends on the relationship between frequencies wq and 
fio. Let tt'< and lu^ denote respectively the smaller and 
the larger of frequencies ujq and Qq- Entering a non- 
dimensional time tn = LOyt into dimensional systems © 
and H14|l we obtain the following value e = w< / uj^ . 

Three cases are possible: 

a) the case of a weak field ojq » CIq, tn — djgi, s — 

b) the case of a strong field ujq « Qq, i„ = flat, e = 

c) the case in which the frequencies are of the same 
order, that is, uiq ^ flo. In this case an additional res- 
onance occurs in the system and the small parameter is 
absent. This situation requires special consideration. 

Clearly, both [E] and [V] are, in general, functions of n 
and m which should be taken into account while carrying 
out estimates. 

11. In problems with initial conditions, the values 
of the coefficients c„|j^q in the expansions of the ini- 
tial functions have order of unity, whilst those which 
do not appear (c„|j^q — 0) are of order of e. As fol- 
lows from Parseval's equality — 1^ coefficient 

c„ must rapidly decreases with the growth of n, thus, the 
first order approximation in contains a finite sum with 
terms having non-zero coefficients c„ of the expansion of 
the initial function. 

It becomes clear from the above that an accurate ac- 



count of all possible internal and external resonances in 
the system, i.e. the analysis of the phase relationships, 
plays a crucial part for obtaining a correct result. For this 
reason, it is natural to refer to this perturbation theory 
as the phase perturbation theory. 

IV. STATIONARY PERTURBATION OF A 
NON-DEGENERATE LEVEL OF THE DISCRETE 
SPECTRUM 

Let us consider a perturbation of a non-degenerate 
level of the discrete spectrum, i.e. the case V = V{f). 
The canonical procedure of averaging is carried out by us- 
ing Hamilton's function (|12|l under the condition ijJ°^„ ^ 
O (e) implying no degeneracy and no close energy levels. 
Simple calculation by means of eq. (|16l) yields 

{Hi}^ 2^ j^c™c„ exp(-?w„„t) 

n,m '^'^ 

fr • \Vkl\^ - _* 

l^k k 

The second approximation to the averaged Hamilton's 
function H12|) is given by 

k 

Acu, = ev,,+e'y2'^-^ (22) 

Both first and second terms on the averaged Hamilton 
function can be rcnormalised by the phase (frequency) 
renormalisation in the original expansion JHJl, i.e. by 
replacing w° by — ui'l + Aujk- In addition to this, 
H^'^^ = 0. This procedure can be performed in any or- 
der of calculations. This means in turn that, instead of 
a standard time interval At ^ 1/e, this approximation 
is valid for exponentially large time intervals which is in 
full agreement with the general theorems of mechanics on 
the behaviour of Hamiltonian systems close to integrable 
systems po). 

The Hamilton function has a diagonal form and the 
coefficients of the quadratic form are corrections to the 
phase (energy) of the unperturbed wave function. 

Equation l|2(J|l of the second approximation for the evo- 

_(2) 

lutionary component has the form 

cfc= -jpr- = -^^^kCk ■ (23) 

"^k 

Then we easily obtain 

= exp{-iAujkt) ■ (24) 
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CoefEcients Ak are determined from the initial condi- 
tions. 

The second approximation to coefficients Ck in expan- 
sion © , obtained by means of formula p9|l , is as follows 



„(2) 



(2) 

Aj^ ' exp(— iAcj/c^) 



eJ2^Al^^exp{-tn^,t), (25) 



mk 



where Vlmk = - ^k,^k ^ ^k + ^^k, and A^,}\ A^^^ 
denote the first and the second approximations to co- 
efficients Afc (it is sufficient to substitute only the first 
approximation to Ak into the second term in eq. H25|l ). 
The second approximation to wave function 5''-^^ is 

(2) 

constructed with the help of coefficients c)j. 



k 



A? - ^^^^ -p( 



x*^exp(-ir!fcO, (26) 



where fifc = uj^ + Aw^, i.e. it is sufficient to restrict 
the consideration by the first correction to the eigenfre- 
quency. 

Provided that the system is in the n ~ th stationary 
state of the discrete spectrum, then A^ — 6kn and we 
obtain from eq. H26|l that 



k 



,0 fc 



exp(— iil„t) 



(27) 



In the case of Cauchy's problem, the system at the 
initial time instant is in a certain stationary state of 
the discrete spectrum, in the s — th state say, that 
is '^{q,t)\i^Q = '^1 and c„|j^q — S^s- Coefficients 
A^'^ , A^'^ are obtained from the relationship 



eq. (|25() . see for example The relationships in eq. 

(|26|l show that it is not correct. This term determines 
the correction to the unperturbed wave function of the 
initial state. The transition probability is determined by 
the second terms in An which is completely absent in 
the standard perturbation theory. This is due to the fact 
that initial condition Cnlt^o ~ ^'"■s is not satisfied. The 
whole coefficient cif' rather than a part of it, as in the 
standard theory, must satisfy this initial condition. In 
the case under consideration, coincidence is occasional 
because the coefficients are independent of time. But 
these coefficients are different in the non-stationary the- 
ory. 

As an example of applying formula (|28|) . we consider 
the problem of the excitation of a charged oscillator by 
an abruptly applied homogeneous electric field e, directed 
along the oscillation axis, ^3|> US- 

In this case it is necessary to solve the problem 



bounded, 



(30) 



where — —ti^A, m denotes the oscillator mass, k is 
the rigidity coefficient, e is the electron charge and V — 
—eex. 

Let us introduce into eq. (|30|l a non-dimensional vari- 
able 1^ — x/a, (a = y/h/nuJ^) and the eigenfrequency 
Lu = y^k/m. The eigenfrequency uj (w = loq) is taken 
as a characteristic frequency loq and the generalised Rabi 
frequency is Hq = eea/h. Assuming the external field to 
be weak, we enter a small parameter e by the relationship 
e = Qq/uj = eea/fku ~ee/ka <C 1. 

An exact solution of the problem of eigenfunctions and 
eigenvalues in terms of the non-dimensional units is given 

by 



E: 



(28) 



From this relationship we obtain An'' — 5ns and = 
5ns — £— TP- The second approximation to the wave func- 

^sn 

tion has the form 



k 

Y^'^-^^leM~^^kt) 



(29) 



^n{x,t) = (2"V7i'a7i!)~i/2gxp 

xi/„ (e - e) exp [-i {n + 1/2 - eV2)] , 

w„ =wo-£V2 = n + l/2-eV2, (31) 



where Hn {z) denotes Hermite polynomials, psj . 

The general solution of problem H3U|) . constructed by 
means of eigenfunctions H31() . has the form 



*(a;,t) = Ecn*n(a;,t) . 



(32) 



Tl = 



This formula is absent in the standard textbooks on 
quantum mechanics. It is important to mention that it is 
adopted in courses on quantum mechanics that the prob- 
ability of transition to this problem is determined by the 
square of the absolute value of the first correction to the 
expansion coefficients c„, that is, by the second term in 



Using the initial condition we obtain the following ex- 
pansion coefficients c„ in eq. 



;"(2"n!)~'/'exp(-eV4), 
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so that the final result is as follows 



where flz, 



ev„ 



This transformation removes 



n=0 

xexp 



xexp [-1 (n + 1/2 - £2/4)] 



(33) 



Performing expansion with accuracy up to terms of 
order of yields 



+e—^^°exp{-it)j exp {-it/2) 

(34) 

where Vl/J^ denotes the eigenfunctions of the unperturbed 
problem. 

Let us now construct the solution of problem (|30l) by 
means of the canonical theory of perturbation. The ma- 
trix elements of the perturbation operator are 



evr. 



+ ((n + l)/2)-i/'5, 



'rn.n+l 



the corrections to the eigenfrequencies are Aujk = — 
the eigenfrequencies of the unperturbed system are ujq = 
1/2, = 3/2, and eq. ^ takes the form 



¥^\x,t) 



+£--^5'?exp(-zi)| exp {-it/2) 



As expected, the solution constructed by the pertur- 
bation theory coincides with the series expansion of the 
exact solution. Let us notice that obtaining solution H34|l 
is simpler than expanding the exact solution in the series. 

The probability of transition to the first excitation 
state is the square of the absolute value of the coefficient 
of *?exp(-i3t/2) and is equal to £"/2. 



V. STATIONARY EXCITATION OF TWO 
CLOSE LEVELS 

In the case of stationary perturbation, it is expedient 
to take into account the principal resonance by renormal- 
isation of the frequency in the original expansion, i.e. to 
present expansion © as follows 



^{q,t) = Vc„(0^°((?)exp(-z^2°i), 



(35) 



the term with svnm from the equations and the effective 
Hamilton function, the latter taking the form 

eHi{c,c*,t) = -ie^ 'i;„™CmC* exp(if]°„i) . (36) 



The presence of two close levels in the unperturbed sys- 
tem, for example with indices a and f3 {uj^ — = eSq) 
means that 17°^ = fl'^ — = e6 {6 — Sq + Vaa — vpjj) 
and leads to the necessity to take into account the de- 
pendence of sum (|36|l on slow time t = et in the process 
of averaging 

eHi = —isVaf}C*^Cf3 exp(?(5r) — ievpaCaC*^ exp(— it^r) 

-ieY, "Vn7nCmC*nexp{-i^nmi)^ (37) 

where notation implies that this sum docs not con- 
tain the diagonal elements (m = n) and the elements 
with m = a,n — P and m ~ j3,n — a. 

The calculation of the first-order approximation yields 
the following averaged Hamilton function 

eHi ~ —ievaf3C^Cf3 e:xjp{iST) 

— ievpaCpCa exp(— i(5r) . (38) 

The Hamilton equations for the evolutionary compo- 
nents Cq, , c/3 have the form 

55i . _ 
e — —levapcp ex-p[idT) , 



C/3 



dill 
'del 



-isvpaCa exp(— i^r) . (39) 



Thus, the averaging procedure "cuts" a two- level sys- 
tem with close levels a and [3 from the whole spectrum. 
It is evident that the case of two-fold degeneracy is a par- 
ticular case of this problem for 5q — Q. The same solution 
has the problem of resonant interaction with an external 
field with frequency w (w — = eSo), the difference 
being only in the resonance type. In the problem of the 
close levels there is an internal resonance. 

The solution of problem (|39|l subject to the initial con- 
dition CqIj^q = 5na is as follows 



1 
A 

A 

0, n ^ a,(3 



[r^i exp(ir22T) — ^2 exp(zilir)] 
[exp(— ifZir) — exp(— if22T)] , 



(40) 



6 A 
where VL1.2 = - ± A 

Representing the first equation in eq. l|in|l in the form 



-ca = £^P^!j^ {A + [exp(-zAr) - 1]} , (41) 
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it is easy to find that the probabihty of being in the 
state become unity {wa — 1) after the time interval r* — 
27r/A. Thus, the oscillations between the levels a and /3 
have period T = 27r/A or frequency uj — A. 

Substituting coefhcients (|40l) into expansion (|35|l we 
obtain that the wave function is a superposition of two 
stationary states with 



(42) 



In the case of a degenerate level 6q = and from eq. 
(|42|l we obtain the correction Acu to the frequency of the 
stationary state = = w° 



Aw 



■4|i;, 



Q/3| 



(43) 



Therefore, all three problems, namely the problems on 
close levels, two-fold degenerate level and resonant in- 
teraction with an external single-frequency field, are all 
solved in the framework of the same approach and yield 
the results coinciding with the traditional one with first 
order accuracy. 

While solving the problem, we determine the condi- 
tions under which the quantum system with a discrete 
spectrum can be modelled, in the first approximation, 
by a two-level system. The main point of this procedure 
is the possibility of averaging Hamilton's function (jSHjl- 
The condition of weakness of the external field is needed 
for this. Then, the non-trivial initial conditions are re- 
quired, at least for one of the coefficients Cq,, c^, otherwise 
the solutions of the homogeneous equations in (|39|l are 
trivial. 

Nowadays, the procedure of solving these problems is 
performed backwards. A two- level systems is first taken, 
then a so-called resonant approximation (rotating wave 
approximation), 13j, .16], is applied to it. At this stage it 
is incorrectly assumed that the resonant approximation 
is also applicable in the cases where the perturbation 
theory is invalid, see for example [l^ . 



VI. NON-STATIONARY SCHRODINGER'S 
EQUATION AS A HAMILTONIAN SYSTEM 

The situation studied in the non-stationary perturba- 
tion theory occurs when operator H can be represented 
as sum H ^ H„ + eV{q,t) (0 < e << 1) of two self- 
adjoint operators. In the adiabatic approximation the 
perturbation operator V{q, t) is not small and depends 
on slow time t = et such that V{q, t) = V{q, r). The so- 
lution should be constructed within the asymptotically 
large time interval t 1/e when change in the pertur- 
bation operator is large. In this case splitting the total 
Schrodinger's operator H into two operators, namely the 
generating (unperturbed) operator and a perturbation 
operator makes no sense. In order to embrace both possi- 
bilities we consider problem Q with the time-dependent 
Schrodinger operator H — H(q,t). 



Let us assume that the stationary problem correspond- 
ing to Q is solvable for a parametric dependence of 
Schrodinger's operator on time and has a discrete spec- 
trum. This means that the eigenfunctions and the eigen- 
values of the problem are given by 



(44) 



with time t being fixed. The eigenfunctions are assumed 
to be orthonormalised as follows 



'lpmiq,t)'lpniq,t)dq = Sr, 



(45) 



where a bar denotes the complex conjugate. 

The existing approximations of Born-Fock and 
Landau-Dykhne ^3 1 suggest that the eigenfunctions 
can be chosen as being real- valued (i.e. no magnetic field 
is assumed) which essentially reduces the applicability of 
the method. In the present study this assumption is not 
needed. 

In the case of weak fields the results of the adiabatic 
approximation of Landau-Dykhne do not coincide with 
the results of perturbation theory |l6j . The approxima- 
tion of Born-Fock is actually not an approximation at 
all since all higher approximations turn out to be of the 
order of the first approximation In addition to this, 
both approximations yield an incorrect factor in front of 
the exponential function 16]. 

Let us look for the solution of the exact problem in the 
form 

^-(5,0 =^c„(i)V'n((?,i)exp|-iy" ri„(z)dz|, (46) 
where 



Unit) = OJnit) + Vnnit), 



LOnit) = E„{t), v„„ = -i I i^n-Q^dq. 



The meaning of this choice of the phase becomes clear 
in what follows. 

Inserting eq. 146|l into eq. Q yields the following 
equation for the expansion coefficients Cm{t) 

Cm{t) = I ^mn{z)dz I , (47) 



where a dot denotes a total time derivative and a prime 
at the summation sign denotes the absence of a diagonal 
components with m = n. The matrix of coefficients Vmn 
has the form 

oo 

Vmn = ~* / '^'rn{q,t) ^^"'^'l'^^ dq (48) 
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and is Hermitian, i.e. Vmn = Vnm- 

The choice of phase indicated in eq. 146|l ensures that 
the sum has no diagonal component which is responsible 
for the principal resonance. If eq. (|46|1 had this diagonal 
component, the sum H47|) would have a small resonant 
denominator. 

Indeed, differentiating eq. H44() with respect to time 
and taking into account that Schrodinger's operator is 
self-adjoint, we obtain 




--— , m ^ n 



(49) 



It is clear that in the case of the real-valued eigenfunc- 
tions -i/jn = Tpm that is, the diagonal elements Vnn = 0. It 
is this fact that is the reason for the real- valued normal- 
isation in the Born-Fock approximation. 

Similar actions in the case when m = n leads to 



the relationship 



dH 



and do not deter- 



dt / dt 

\ / nn 

mine the diagonal matrix elements. In the case in which 
Schrodinger's operator depends on time r in terms of 
the set of functions ^,;(t) {i — 1, 2..N) , elements w„„ de- 
termine the topological adiabatic Berry phase whose 
value does not depend on the evolution time and is deter- 
mined only by a closed contour in the parameter space. 

Equation H49I) indicates three cases allowing the de- 
velopment of the perturbation theory. In the adiabatic 

~ , , dH ( dH\ fdC 

case H = H{^{t)), so that — = ^ I I ( 

In the case of the non-stationary perturbation theory, 
Schrodinger's operator has the form H — Hq + £V{q,t)- 
Finally, in the case of the adiabatic perturbation theory 
H ^ Ho + eV{q,T). 

The system of equations (I47|l , along with the complex 
conjugate one, is Hamiltonian (in the classical sense) hav- 
ing the following Hamilton function 

H{c, C*,t) = -i^'v„in{t)Cn 



X exp {-ij ^mn{z)dz } (50) 

which describes the classical distributed system. The ma- 
trix of coefficients Vmn is Hermitian which ensures that 
this system is Hamiltonian and in turn enables one to 
apply the phase perturbation theory. 



function H50() can be cast in the form 
£i?i(c, c*,t, r) ^ -ie^'vmniT)c„ 



XC*^exp{iJ QrnniT)dt } , (51) 

where t = e and = — ^n- 

The canonical form allows us to convert the evolution- 
ary equations by means of the formulae 

(2) (c, c* , r) = eili (c, c* , r) + e^H2 (c, c* , r) , 

..H..,ft^-((t)(^l^)>, 

where -ff*^^-* denotes the second approximation to the av- 
eraged Hamilton's function, whilst c — (ci,C2,..) and 
c* = (c*,C2, ..) are the evolutionary components of vari- 
ables c and c*. 

Averaging expression l|51|l along the generating solu- 
tion (cfe = const, T = const) , we obtain Hi = (Hi) = 

_• 

which in turn implies that Cfc = 0. The latter result is 
the adiabatic theorem of Kato which is obtained in 
fact without calculations, cf. ((30J) for the proof. In the 
classical sense, the evolutionary components Ck of the 
orig inal variables Ck are the adiabatic invariants, see j^, 
[20j, i.e. they retain the initial values for the asymptotic 
time interval t ~ l/e. For deriving this result it is neces- 
sary to assume that Qmnir) ^ 0{e), i.e. the system has 
no degeneracy, there are no close levels and the levels do 
not intersect during the evolution time. 

In the adiabatic (first) approximation, the solution of 
Schrodinger's equation has the form 



vI/W(g,i) = ^c«V«('Z,i) 



n=0 



X exp i ^^J ^mniz)dz ^ . (53) 

Under rather general assumptions, Los's theorem 
renders estimates for the difference | '^{q,t) — ^^^\q)\ < 
Ce, where C is a constant independent of e for time in- 
terval t l/e. 



VIII. POST-ADIABATIC APPROXIMATION 



VII. ADIABATIC APPROXIMATION 

We assume that Schrodinger's operator has the form 
H = Hq + V{q, ^{t)), where t — et denotes slow time. 
Then the matrix elements Vmn ~ ^^id Hamilton's 



In order to construct the second (post-adiabatic) ap- 
proximation we make use of relationships in eq. H52|l . 
Simple calculation yields 



H, 



k 



\Vkl\_ 



(54) 
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so that the second approximation to the averaged Hamil- 
ton's function H^^^^ has the form 

= -^£2^A^]fe(T)cfcC^ (55) 

A; 

Hamihonian equations with Hamihon's function H55|l 
for the evolutionary components Ck are integrated easily, 
to give 

Ck = Afcexp |-ie^y Af2fc(z)dz| 

= Ak exp (-iafc) . (56) 

Integration constants A^, are determined by means of the 
initial conditions. 



Let us notice that the phase of coefficients Ck could be 
included into the original expansion gnil, then we would 
obtain H2 = 0, i?^^) ^ 0. 

The second approximation to the expansion coeffi- 
cients in eq. (|46(l is constructed by means of formulae 

m 

(2) 4 (2) r ■ \ 
4 - ^fc exp(-iafe) 

_ ^y-^ l^^kvn_j^^ exp(-iQ!m -I- iVLkrat), (57) 

With the help of coefficients c). we obtain the second 
approximation 5'(^)(g,t) to the solution of Schrodinger's 
equation 



Afeexp(-iafe) - gy^^T^^^m exp( 



km 



exp 



—i J ilk{z)dz 





(58) 



It is necessary to mention that in eq. (|57(l we can limit 
our consideration to the terms in the sum by the first 
approximation A^n with respect to e. 

When Cauchy's problem is studied, the system at the 
initial time instant is at a certain stationary state, say s — 
th, of the discrete spectrum of the unperturbed problem 
with Schrodinger's operator Hq, i.e. ^'('?,i)lt^o ~ 
In contrast to the stationary case we can not take that 
Cnlt=o ~ '^"S' si'^ce the expansion is carried out in terms 
of the eigenfunctions of the perturbed problem ipniq, t) — 
'0n(9jC(O); where ^(t) denotes parameters determining 
the dependence of the perturbation on time. With this 
in view, we additionally assume that the switch fulfills the 
conditions ^(0) = ■^(0) = 0. The problem can be solved 
under other conditions which implies an instantaneous 
switching of the perturbation followed by its adiabatic 
change. 

Therefore, we take that the equation for coefficients Ak 
has the following form 



Ck 



t=o 



A, 



Vkmir) 



^kmir) 



Ari 



(59) 



Obviously the matrix elements Vkm ~ ^(''')i thus under 
the adopted conditions we obtain Aks — Sks, that is 



„(2) 

^ks 



Vk. 



Sks exp{-iak) - e—— exp{-ias + iQkst) ■ (60) 



n 



ks 



Finally, the second (post-adiabatic) approximation '^^^^ 



for the wave function is as follows 
¥^^q,t) = J2[AkcM-iak) 



, Vkm 
^km 



Amexp{~iam - i^kmt) 



(61) 



X ^/kcxp 



-i J flk{z)dz 





It can be proved, see that the estimate 

\'${q,t) — 5'(^)(g,t)| < Be-^ is valid for the time interval 
t ~ 1/e. Papers by physicists do not take into account 
the boundedness of the time interval in which an approx- 
imate solution approximates the exact solution. As one 
can see from the forthcoming examples, this interval, in 
general, can not be enlarged. 

The problem of constructing the adiabatic approxima- 
tion and taking account of the transition through the 
virtual levels was posed by Dykhne: "In order to obtain 
the correct factor of the exponential function it would 
be necessary to take into account all higher approxima- 
tions of the perturbation theory, all yielding results of 
the same order. In practice, this is, of course, not feasi- 
ble. The obtained formulae give answers to the question 
of calculating the probability of transition of a quantum 
system to an "adjacent" level. As for transitions to more 
remote levels, then the transitions through virtual levels 
may compete with the considered process of the " direct" 
transition. However, this question needs an additional 
investigation", Nevertheless, this problem has not 
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been solved so far by the existing methods of perturba- 
tion theory. 



IX. QUANTUM LINEAR OSCILLATOR IN A 
VARIABLE HOMOGENEOUS FIELD 



In order to compare the obtained results with the 
known ones, let us study the following problem having an 
exact solution. The situation considered is the motion of 
a particle in the field of a parabolic potential subjected 
to a variable external force, i.e. 



ih 



dt 



-1/4^1/2, 



9 2 

muj q 



eeaf{yt)q 



■^{q.t) 



*(g,OI*=o='^^^ « ^ exp 
*(g,Olg^±oo - bounded. 



(62) 



where m and e denote the mass and charge of the oscil- 
lator, eo is amplitude of the electric field, denotes a 

1 /2 

characteristic time constant of the field, a — {h/muj) ' 
denotes a characteristic length scale. Clearly, the initial 
state is the main state of the free harmonic oscillator. Let 
us introduce the non-dimensional time tn = cot, the non- 
dimensional coordinate x = q/a and the non-dimensional 
force amplitude Si = eeoa/hio ~ fl/uj. Problem (|62|l in 
the terms of non-dimensional variables takes the form 



*(x,t). 



*(a;,i)lt=o-'^"'^'exp(-a;V2), 
*(a;, t)\^^+oo ~ bounded , 



(63) 



where e = i^/uj, ^{et) = eif{et) and the non- 
dimensional time is denoted by t. 

The exact solution of this problem is given by, cf . [s^l , 



In order to construct expansion (|64() in the adiabatic 
case (0 < £ << 1) we integrate eq. (I65II by parts three 
times and take into account that ^(0) = ^(0) = <^(0) = 0. 
The result is 



i + ii-i 



exp(ji) ~ 2 I exp(iz) ^ {z)dz . 



Taking M = max 
the following estimate 



(66) 

^(t) in time interval [0,T] we obtain 



exp(iz) ^ {z)dz 



^ J \^{z)\dz MT ^ e-'T . (67) 




It is clear from this equation that one can neglect the 
latter term in eq. (|56|l within time interval At ~ T ~ 
l/£. For this reason, it is necessary to keep the second 
order term ^ in eq. H66|l for evaluating integral (|64|) and 
omit it by substituting without integration. 

With the same accuracy we evaluate the following in- 
tegral 



e 



e 



e 



i J dzexp(~2iz)S (z) = i J -^dz — - — + i J ~^dz 
00 

(68) 

and the expansion of the exact solution H64|) as a series 
in time interval t ^ l/£ 



it 

expi--- 




O e' 



(69) 



Next, we construct the expansion of the exact solution 
(|64|l in the case of the harmonic non-resonant perturba- 
tion {£, {t) = £1 sinvt, < £1 << 1,1^^0 (£1) , V ^ 1). 
In this case 



^o{x,t) = TT cxp < + dz5^(z) exp(-2iz)- 
I 

V2x6{t)e^p{~it)'> , (64) 



S{t) ^ ——^ J ^{ez) ex.p{iz)dz (65) 


for arbitrary values of the parameters e and £1 . 



S{t) = -£j 



V2{iy^ - 1) 



[{i sin vt — V cos vt) exp(it) + v] 



(70) 



I J dzS^iz) exp(-2it) = -» 4(^2^_ + O i^i) C^l) 


The terms omitted in eq. (|71f) are uniformly bound 
such that approximation H71|) is valid at any time instant. 
Inserting eqs. H70() and eq. H71() into eq. (|64|l yields the 
following expansion of the exact solution 
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*o(a;,t) 



-1/4 



exp ■ 



it 
2 



^1 



^4(1/2-1) 



+ x- 



[isuivt— cos J't) + exp(- 



^'-1 ..2 1 



5*1 (i sin vt — V cos I'f) 



it 

exp 1 -- 



leiv 



V2i 



, .3t 
^-vl/iexp^-.- 



exp 



el 



where -00 and i/'i are the eigenfunctions 
of the unperturbed Schrodinger's operator 
{^{t) — El cost, < £i << 1). In this case 



6{t) 



lei 



(exp (it) sint + t) 



I y (5 (z) exp(--2iz)dz 




(72) 
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1 / sin2t 1 9 
2 



t^ exp(— it) sint 



so that within non-dimensional time interval t ~ 1/e the 
expansion of the exact solution with accuracy up to the 
first order of smallness has the form 



4-0(2;, t) 



-1/4 
TT ' exp 



it 
2" 



Y 



Eit 



X. 



ia;^exp(-it) + ) cxp(-2it) |>(73) 



CHARGED LINEAR OSCILLATOR IN AN 
ADIABATIC HOMOGENEOUS FIELD 



In order to demonstrate application of the suggested 
theory, let us solve the problem of motion of a particle 
in the field of a parabolic potential subjected to a vari- 
able external force. Let e = v/co he a small parameter 
(0 < e << 1) in the adiabatic case. Let us also take that 
parameter Si is of order of unity. The perturbation op- 
erator ^{£t)x can not be taken to be small in any way, 
cf. 0, that is, the problem of the perturbation theory 
can not principally be solved by the methods of spectral 
analysis of operators. 

To make calculations more transparent, it is worth- 
while carrying out the calculations from the very begin- 
ning rather than to use the resulting formula l|t)l|) . For a 
fixed time instant t problem (|44ll is written in the form 



1^2 

2di2 



*„(a;,t) = En{t)-^n{x,t) 



(74) 



and has the following solution 



^n{x,t) =(2"V^n!) 



n-1/2 



exp 



-{x-aet)r/2 



xHn {x - ^(et)) , 
w„ = n + 1/2-^2/2^ 



(75) 



where Hn{z) denote Hermite polynomials [13, |28l |. 

Using the recurrent relationships for Hermite polyno- 
mials it is easy to calculate the matrix elements Vmn{t) 
for m ^ n 

Vmn{t) = - (iif^j [Vn+ l(5„i,„+l - VnSm^n-i] (76) 

where a dot denotes the total derivative with respect 
to time, so that ^ ~ e and Sm,n denotes Kronecker's 
delta. Since the eigenfunctions l(7S|) are real-valued, the 
diagonal matrix elements Vnn = 0, i.e. Sl„ = ujn = 
n-f 1/2-^2/2. 

With the help of eq. ifTHjl we can easily calculate values 
Arife = Aojk = — ^/2 and coefficients Ck, see eq. (IKI)|> 



Cfe = ^fcexp<i/ i'^i£z)/2 dz 



(77) 



In the case under consideration (Hi) = and Hi = Hi, 
d{Hi} 



then we find {Hi} and 



del 



{Hi} = -i^'^^^^^^^c„cj;,exp[ia;,„„(T)t] 



d{Hi} 
del 



,Vkn{t) 



:„exp[iWfc„(T)t] . (78) 



tOkn (t) 

Making use of relationships (|78|l we obtain 



„(2) _ 4(2) 



A):' exp 



■ t 

i / — dz 
J 2 

. 



E 



, Vknjt) 



An exp 



I / —dz + iUknt 




(79) 



13 



The equations for determining the integration con- 
stants have the form 



,(2) 



= Sko = Ak — ^ ' 



, f Vknjt) 



An 



(80) 



Assmiiing, as above, that ^(0) ~ Vkn{Q) — 0, we obtain 
Ak = 5ko and 



(2) X 

Ck = Ofco exp 



" * '2 

i f —dz 


VkO 


exp 


J 2 


Wfco 


. 





dz + iukot 



As follows from eq. (|7S|l Vko = —i£./V^, thus 



„(2) 



v2 



exp 



(81) 



(82) 



Inserting these equalities into expansion (|46|l yields the 
post-adiabatic approximation 



it 

X exp<; -- 



V2 



vI'i(x-0 



—dz + i f —dz 
2 ./ 2 



(83) 



XI. ADIABATIC PERTURBATION THEORY 

In the case of the adiabatic perturbation theory, 
Schrodinger's operator is as follows H = Hq + eV{q, r) 
which allows us to carry out calculations by using the 
results of Sec. 6. 6. However, as mentioned in Sec. 6.2 it is 
more convenient to apply the formalism of the stationary 
phase perturbation theory. 

Indeed, by casting problem H62|l in terms of the eigen- 
functions of unperturbed Schrodinger's operator Hq we 
find the effective Hamilton function 

eHi {c,C*,t,T) ^ ^^^X! (*^mriO ' 

71, m 

(86) 

where VmniT) denote matrix elements of the perturba- 
tion operator which are calculated by means of the un- 

t 

perturbed eigenfunctions, further — tjJJ+e J VnniTjdr 



and ri°j„ = il^j — ri", denote the eigenvalues of the 
unperturbed problem. 

In this case formulae (|52|) remain valid since the depen- 
dence of the effective Hamilton function on slow time r 
is observed only in terms of the third order of smallness. 

For problem H62|) we have 



which coincides with the expansion of the exact solu- 
tion H69|) within time interval Ai ~ 1/e, approximation 
(|83|l not being applicable for increased times. 

As follows from eq. (|83|l . the probability of the os- 
cillator excitation is equal to ^^/2, i.e. it is zero at ex- 
treme points of function ^(i). For a Gaussian distribu- 
tion, ^{t) ~ exp(— T^) and the excitation probability has 
its only maximum at r = 1/ ^/2. 

In the traditional Born-Fock approximation, [TSl , , 
the expansion coefficient corresponding to cq. ()79|l is 
given by 

t 

Cfcn = 4n+ / — AtTV f exp ij UJkn{t")dt" dt' , 





UJknit') \dt J 



where l^(x,t) = —£^{t)x. Carrying out simple manipula- 
tions we obtain 



CfeQ = ^fc( 



fcO 



4i 
V2 



iexp{it')dt' . 



(85) 



*51(x) = (2"^/^n!)-i/2exp [~xy2] i/„ (x) 
o;^ = n + 1/2, = =n+l/2 
ei^ir) \ ln + 1 



'^mn ('^) 



VnniT) = 
1 

2 



Aujk 



V2 
e\e{r) 



+ 1 



(87) 

with parameter ei being a small value of order of e. 
By virtue of the latter relationship in eq. H87|l we eas- 
ily find the second approximation to the expansion coef- 

(2) 

ficients c„ (under the conditions ^(0) = ^(0) = 0) 



exp 



IS' 



' ' ^dz 



3fe0 



and the second approximation to the wave function 



Comparing the latter equation with eq. H82|l indicates 



that the factor exp 



^/ (e'W/2) 



dz 



is absent in eq. 



(|85|l . Equations H85|l and (|82|l coincide at time instant 
t ~ 1 , when a change in f {et) in the integrand can be 
neglected and there are no conditions ^(0) = ^(0) = 
under which approximation l|82|) is valid. 



exp 



-- — \- i [ — dz 
2 J 2 





(89) 



which coincides with adiabatic expansion (|83|l when we 
take into account a small factor at ^ due to parameter 

£l. 
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XII. HARMONIC EXCITATION OF A 
CHARGED OSCILLATOR. NON-RESONANT 
CASE 



In this case Hi = (Hi) = 0, so that Hi = Hi- The cor- 
rection of the second order H2 to the averaged Hamilton 
function is as foUows 



To demonstrate the way of constructing solutions in 
the case of non-stationary perturbation theory, we con- 
sider the case of a harmonic external field which is fre- 
quently encountered in practical applications. In this 
case function £^{t) in eq. H62() takes the form ^(t) — 
eisinj/i and we can adopt that v 1 and 7^ 0{e). 
In other words, we consider the non-resonant case and 
parameter Si is taken as having the order of smallness of 
parameter e. 

Let us carry out the corresponding calculations in three 
ways: first, stationary phase perturbation theory, second, 
non-stationary phase perturbation theory developed here 
and finally, traditional method of p^ . 

In the first case the matrix elements of the perturba- 
tion operator, calculated by means of the unperturbed 
eigenfunctions, have the form 

Vninit) = £i(wm„ exp(wt) - exp(-wi)) , 

^mn = [Vn- + I5m,n+1 + V^^m.n-l] ■ (90) 

2v 2 

The effective Hamilton function is given by 



Ho = -i 



4(j/2 - 1) 



(92) 



Next we find the second approximation for variables c„ 



(2) 



V2(l - ly^) 



— %—f= [isini/i — vco'&vi\ exp(it) . (93) 
v2 



£liJl(c, C*, i) = '""m" *5^P + ^)^] '^^'^ 



iY] 'v^n exp [iiuj'^n ~ t^yi CnCl, . (91) 



n,m 



Taking into account these relationships we obtain the 
second approximation ^''^^^(x,/:) to the solution of the 
non-stationary problem 162|l 



lei 



\/2(i 



-'^ i{is\a.vt — vcosvt) 



exp 



-^_^v,,exp(^-^-,^exp 



4(1/2 _ 1) 



(94) 



This result suggests that the spectrum remains equidis- 
tant and is only subjected to a common shift Alu — 
el/ A {v^ — l)- In this case Aw < for i/ < 1 (cj < wq, 
Lo being the frequency of excitation force) , and Aw > 
for 1/ > 1 (w > Wo). The shift Aw at w ^ corresponds 
to effective elevation of the bottom of the potential well. 
In the case of a high frequency external field (v >> 1) 
Aw — el/Av"^ that coincides with the effective potential 
energy. Hence, the constructed expansion l|94() coincides 
with the expansion of the exact solution (|IX|I . 

Let us now construct a solution by using formulae of 
the non-stationary phase perturbation theory. In this 
case the matrix elements (|48|l calculated by means of 
eigenfunctions (|75() of the instantaneous Schrodingcr's 
operator have the form 

Vmnir) = -= {VnTTSjn,n+l " VnSm,n-l) (95) 



or 

Vmn{t) = uj^„exp(wi) -|-u^„exp(-ii/t), 

f^mn = ~*7r7^ {Vri~\^Srn,n+l ~ VnSra.n-l) (96) 

2\/2 

As the eigenfunctions are real, then u„„ = and ri„ — 

^1 e 

The effective Hamilton function is as follows 

OG 

eiHi{c,c*,t) ^ -iei ^ '{v^^exp [i{Lu^„ + v)t] 

m,n— 

+v'Ln exp [t(t^™„ - i^)i] jCnC ■ 

(97) 

It is evident that Hi = (Hi) = 0, so that Hi = Hi. 
The correction H2 to the averaged Hamilton function, 
calculated by formulae H52|) with the help of eq. (|96|) . is 
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given by 



'4( 



where A^^^^ and A^^'^ are the first and the second approx- 
imations to coefficients with respect to parameter £i. 
(98) As the first approximation we can take A^.^^ = 5ko, then 



The second approximation for the evolutionary com- 
ponents of variables is found from Hamilton's equa- 
tions with Hamilton's function H^"^^ — eiHi + e\H2 and 
has the form 



Cfc = Ak exp 



4(z/2 - 1) 



Constants Ak should be obtained from the initial condi- 
tions. 

(2) 

The second approximation cj. to the original variables 
Cfc is as follows 



Ak-eiY, vl^A^ 

n 

exp [i {iOkn " t] 



exp [i {uJkn + v) t] 



LUkn + V 



UJkn — V 



exp 



(100) 

For determining the integration constants Ak from eq. 
(|100|l we obtain the following equation 



A 



(2) 



5ko - i 



^/2(l - 



-Ski 



(102) 



where it is taken into account that w^q = —ieiv5kil2\/2 



(99) and the final expression for Cj. takes the form 



eiv £ii'exp(it) 
-Okl + 1- 



V2(i/2-l) 



<Ski [cos vt — iu sin ut] } exp 



4(z/2 - 1) 



(103) 

Then we find the expression for phases in expansion l|46() 



r ^2 I 

-i I LOn{t)dt = i—t - i{n + -)t + O(e^). (104) 



c/clt=o — Sko — A 



(2) 



£1 



(iVkn + I') ^ + {uJkn - v) ^ 



(101) 



Inserting eqs. p03|l and p04|) into expansion (|46|l we 
finally obtain 



*(')(a;,i) = 



- 1) 



1 (i sin vt — v cos vt) 



it 



isiv ( 3i 

X exp I — — I H — Wicxp — i — 



(105) 



exp 



^4(z/2- 1) 



r 



coinciding with eqs. 1)94(1 and IjlXI) . 

Finally, taking into account that the second approx- 
imation of the present analysis coincides with the first 
approximation of the traditional approach, we obtain, 
for the case ak = Ck, by means of the standard formulae 
of the non-stationary perturbation theory, that 



a-k 



JO) 



,(1) 



Sko + i 



£1(5, 



kl 



- 1) 



(i sinj/t — :/cos i^i) exp(it) . (106) 



Comparison of eqs. p03l) and (|106|l allows us to indicate 
a number of inaccuracies in the standard courses. Firstly, 
a phase multiplier eyi\)\—i£\v^t/ Aiy^ — 1)], which is of 
crucial importance for investigation of the coherent pro- 
cesses, is absent in expression (|106f) . It can be neglected 
within a non-dimensional time interval i ~ 1, but not 



within asymptotical intervals t ^ \/e. Thus, approxima- 
tion 1106|l and in turn the whole solution is valid only 
within this small time interval. 



Secondly, only the first approximation a]^^ = 6, 



kO 



rather than ak is subject to the initial condition (|101|l . 
This explains the absence of the term i£li//\/2(j/2 - 1). 



,(0) 



Bia^ko ^^'^ subjecting 



Indeed, by assuming aj. - — Oko 

the whole coefficient ak to the initial condition afel^^Q 
Sko we find 



-(0) 



V2(;^2 _ 1) 



Ski 



(107) 



Finally, as follows from solution (|105|) or l|IXp it is this 
absent correction (rather than a[,^^) that determines the 
probability of transition to the excited state. 
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XIII. HARMONIC EXCITATION OF AN 
OSCILLATOR. TRANSITION THROUGH A 
RESONANCE 

Let us consider excitation of an oscillator by a weak 
resonant harmonic field V {x,t) ~ —eix cosvt. In this 



case 1 — i^ = e (0^e<l) and parameter ei is a small 
value. 

The effective Hamilton function constructed by means 
of the eigenfunctions of the unperturbed Schrodinger's 
operator has the form 



eHi{c, c*,t) — i-^^ ^ { \Aic„c*_i exp [-iet] + \/n + lc„c*^i exp [ist]} 



+ i^-^ ^2 {Vnc„c'^_i exp [—i {v + 1) t] + \/n + lc„c*+i exp [i [v + l)t\] . (108) 



r 



In the first approximation we obtain the averaged and solution 'I' (x, t) in the first approximation 
Hamilton function 



sHi = i — ^ "S^ {\/nCn c*_i exp [—iet] 
2v2 V 



Cn c* +2 6xp [iet] } 



(109) 



k=0 



exp 



expi-^) (2^-V^fc!)-V2 



Hk (x) exp 



-i[k+-]t 



(11 



and the equation of first approximation for the evolution- 

ary components Cfc of the original variables takes the Carrying out summation in eq. ^ by means of the 

generating function for Hermite polynomials [2§| 

oo ^ 

exp{2xz- z^) ^Y^^Hkix), (113) 



form 



Ck= « 



2\/2 



Vfc + lcfe+1 exp {—iet) + Vkck-i exp (iei) 



(110) 

It is easy to prove by direct differentiation that the 
solution of the equation 



k=0 



we obtain the final solution for the case of exact resonance 



Cfc = i- 



V2 

has the following form 



Vk + Icfe+i exp {—iet) + Vk exp (iet) c^-i 



'^{x, t) = TT 'i exp ■ 



it x"^ 
2" ~ Y 



eit 



e\t . , I e\i \ / r\ \ 
IX — exp(— it) + — I cxp(— 2ii) 



co{t) = exp l^^J dzS'^{z) exp{—2iez) + (t) exp{—2iet) > , 



(114) 



Ck{t) 

6{t) 



This solution coincides with the expansion of the exact 
solution (|7S1) . The probability of excitation of the oscil- 
lator has the form of a Poisson distribution 



■co(i) , 



V2 ^(^) '5^P(*^^)^^ ■ 



(111) 



(n)" 

Wn(t) = |c„(i)|^ = — |-exp(-ri). 



(115) 



Using this solution one can investigate the transition 
of the system through the resonance. 

We restrict our further investigation to the case of ex- 
act resonance (e = 0). Using relationship Hlll|) we find 
coefficients 



where n = {eit/2^/2)' 



Ck{t) 



1 



/ ieit\ 
\2V2) 



exp 



eit 
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